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Cosmological parameters constraints from galaxy cluster mass function 
measurements in combination with other cosmological data 
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Abstract — We present the cosmological parameters constraints obtained from the combination of galaxy cluster mass 
function measurements (Vikhlinin et al., 2009a,b) with new cosmological data obtained during last three years: updated 
measurements of cosmic microwave background anisotropy with Wilkinson Microwave Anisotropy Probe (WMAP) 
observatory, and at smaller angular scales with South Pole Telescope (SPT), new Hubble constant measurements, baryon 
acoustic oscillations and supernovae Type la observations. 

New constraints on total neutrino mass and effective number of neutrino species are obtained. In models with free 
number of massive neutrinos the constraints on these parameters are notably less strong, and all considered cosmo- 
logical data are consistent with non-zero total neutrino mass Era„ ~ 0.4 eV and larger than standard effective number 
of neutrino species, Neff ~ 4. These constraints are compared to the results of neutrino oscillations searches at short 
baselines. 

The updated dark energy equation of state parameters constraints are presented. We show that taking in account sys- 
tematic uncertainties, current cluster mass function data provide similarly powerful constraints on dark energy equation 
of state, as compared to the constraints from supernovae Type la observations. 
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INTRODUCTION 

The measurements of galaxy cluster mass function 
give one of the most sensitive method to measure the 



dark energy equation of state 


fe.g., IStarobinskvl 


1998 


Haiman et al. , 2001 ; Weller et al. 


,20021 Wans et al.. 


2004 


Maiumdar, Mohr, 2004). Current observational data 



provide measurements of cluster mass function accurate 
enough to obtain powerful c onstraints on the parameters 



of cosmoloaical model (e.g., Borgani et al. 2001t Henrv 


2000 


, 20041 Reiprich, Bohringer 2002; Vikhlinin et al. 


2003 


: Schueckeret al., 2003t Voevodkin, Vikhlinin 


2004 


: Mantz et al. 


20081 2010al Vikhlinin et al. . 2009b 


Vanderlinde et al., 


20101: Vikhlinin. 2010h. 



In this work the confirmation of the existence of the dark 
energy was, for the first time, obtained using new indepen- 
dent method, based on the measurements of the large scale 
structure growth rate, not only on the measurements of the 
geometry of Universe. The joint analysis with other cos- 
mological data allowed to significantly improve the mea- 
surements of the dark energy equation of state parameters 
and also to improve constraints on other cosmological pa- 
rameters, e.g. total mass of light neutrinos. 

In the last years a significant amount of new cosmo- 
logical data were published. These inclu de seven-year 



data of CMB observations with WMAP (ILarson et al 
201 lb. and CMB observations at smaller angular scales 



The strongest to date constraints on cosmological pa- 
rameters using the data on meas urements of galaxy clus - 
ter mass function were obtained in Vikhlinin et al. I d2009bl) . 
These constraints appears to be similarly powerful as 
compared to the constraints from most recent cosmic 
microwave background (CMB) measurements made by 
WMAP, observations of baryon acoustic oscillations and 
supernovae Type la. These new constraints are independent 
and have different cosmological parameter degeneracies. 



dReichardt etafl 120091: iBrown et all 120091: bunklev et al 
2010: lKeisler et all 1201 II) . Using the calibration of super- 



novae la absolute luminosities the mea surement of Hub- 
ble co nstant was significantly improved (Rie ss et al. i l2009l 



2011). Also the new data on observations of s uper- 
novae la ( Hicken et al. , 200St Amanullah et all 2010l) . and 
baryon acoustic oscillations jPercival et al. , 2010l) were ob- 
tained. Cosmological constraints from the joint analy- 
sis of_these_andothe£sc presented 



Komatsu et al.l (1201 ll) . iKeisler et alJ (120111) and other 



e-mail: rodion@hea.iki.rssi.ru 



works. 

In this paper we present the cosmological parameters 
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constraints, obtained in result of joint an alysis of galaxy 
cluster mass function measurements dVikhlinin et al 



2009albT) and re cent cosmological data discussed above. 
As compared to Vikhli nin et al. I d2009bl) . wider set of cos- 
mological parameters is considered. To calculate joint 
likelihood functions in cosmological parameter space we 
used Markov Chain Monte Carlo technique (see, e.g., 
Lewis. Bridlel 120021) . In order to include galaxy cluster 
mass function cosmological data in these calculations, cor- 
respondent software was developed, which we provide for 
public use. 

In our work the new significant constraints on total neu- 
trino mass and effective number of neutrino species are pre- 
sented. In order to test the possibility of the existence of 
light sterile neutrinos with masses near 1 e V, which were 
sugge sted to explain t he results of LSDN dAguilar et all 
200lh and MiniBooNE dAguilar- Arevalo A. A. et al.Ll2010l) 
experiments, and also recently discovered so-called reac 



tor ne utrino anomaly dMueller et all 1201 It iMention et al 



201 ll) . we consider the models with free number of mas- 
sive neutrinos. Also in this paper we give the updated con- 
straints on the parameters of the dark energy equation of 
state. We show, that taking in account systematic uncertain- 
ties, dark energy equation of state constraints based on the 
existing cluster mass function data and supernovae Type la 
observations are comparably powerful. 

COSMOLOGICAL DATA 
In our work the data on cluster mass function mea- 



surements are taken with no changes from IVikhlinin et al 
(2009b). In this work cluster mass function was measured 
using a sample of 86 massive galaxy clusters with masses 
measured usin g Chandra X-ray obser vations with about 



10% accuracy dVikhlinin et all l2009ak . Distant (z « 0.4 



0.9) clusters in this sample were selected in 400 square de- 
gree X-ray g alaxy cluster survey, based on ROSAT PSPC 



pointed data dBurenin et all 120071) . Clusters in local Uni 



verse (z < 0.2 ) were selected from RO SAT All Sky Survey 



(see details in [Vikhlinin et alll2009al) . Confidence regions 



for different cosmological parameters were obtained from 
the analysis of likelihood function, which was calculated 
on a grid of cosmological parameters. The results of these 
calculations are available on WWW 1 . This cosmological 
data set is designated below as CL. 

In our work we used also new cosmolog ical data, which 



were pu blished after the issue of paper by IVikhlinin et al 
d2009bk . The most significant improvement was achieved 
in the measurement of Hubble constant, due to the cali- 
bration of supernovae la absolute magnitudes, which was 
made using the Cepheid observations in SNe la host galax- 
ies. After that from the observations of nearby SNe la the 
measuremen t Hq = 73.8 ± 2.4 kms~ 1 Mpc" 1 was obtained 
dRiess et all 1201 11) . These cosmological data are desig- 



nated below as Hq. This is much better measurement a s 
comp ared to Hp = 72 ± 8 kms~ 1 Mpc~ 1 dFreedman et al 



200 lk used by IVikhlinin et ail d2009bk . and also as com- 
pared to the meas urement Hq = 74.2 ± 3 . 6 km s~ 1 Mpc~ 1 
dRiess et alll2009k used bv lKomatsu et al.1 d201 lk . 

In our work we used s even-year da ta of C MB ob- 



servations from WMAP dLarson et all 120111) . For 



these data likelihood were calculated using the soft- 
ware from LAMBDA archive, version 4. 1 2 . The- 
oretical CMB anisotropy spectra were calculated us- 
ing CAMB software, version January 2011, where re- 
combination was c alculated using RECFAST software 
dSeager et al. , 1 19991) . version 1.5, which incorporate a 
number of improvements in hydrogen and helium re- 
combination theory, obtained during the last year s (e.g., 
Dubrovich. Grachev! 120051: IChluba. Sunvaevl l2006i 12007 



2009 



2007 




Sunvaev. Chlubal 120091: iKholupenko et al 
This cosmological data set is designated as 
WMAP7. 

Also we used recently published data on CMB 
anisotropy at smaller angular scales obta ined with South 



Pole Telescope (SPT, iKeisler et al 



201 11) . The data were 
The contribu- 



taken from the web page of the project 3 
tions of "poisson" and "clustered" sources were added to 
the theoretical CMB anisotropy po wer spectra accordin g 
to the directions in §4.1 of paper by IKeisler et alj ( 1201 lk . 
The templates for these components were also taken from 
the project web page. This data set is designated below as 
SPT. 

In addition we used the data on baryon acoustic oscilla- 
tion measurements from the data Sloan Digital Sky Survey, 
Data Release 7, and the data of 2dF GRS survey. These data 
were taken as Gaussian priors for the ratio of sound hori- 
zon at the baryon-drag epoch to effec tive distance measure 
for two redshifts z = .2 and z = 0.35 ( Percival et al. , 2010t 



Komatsu et all 1201 II) . This data set is designated as BAO. 



In order to study the dark energy equation of state we also 
used t he data on SNe la observ ations, the UNION!"" compi- 
lation ( Amanullah et al. , 2010k . This data set is designated 
below as SN. 



LIKELIHOOD CALCULATION 



In Vikhlinin et al 



(2009b) likelihood functions for 
galaxy cluster mass function data were calculated at cos- 
mological parameters grids, since in that case these cal- 
culations could be performed very effectively. In order to 
use these likelihood calculations in Markov Chains Monte 
Carlo simulatio ns in multidimension al cosmological pa- 
rameters space dLewis. Bridlel 120021) . we wrote a mod- 



ule for CosmoMC 5 software, where likelihood in arbitrary 



http://hea.iki.rssi.ru/400d/cosm/ 



2 http://lambda. gsfc.nasa.gov/product/map/dr4/likelihood_get.cfm 
3 http://pole. uchicago.edu/public/data/keislerl 1/ 
4 http://supernova.lbl. gov/Union/ 
5 http://cosmologist.info/cosmomc/ 
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point of parameters space is calculated using simple inter- 
polation on parameters grids. In spite of ideological sim- 
plicity, the realization turn out be somewhat cumbersome 
due to many dimension grids. Therefore, we think it may 
be useful to provide this module for public use. This soft- 
ware can be found at web page of CCCP project 6 . 

Using this module one can easi ly include likelihood for 
galaxy cluster mass function from lVikhlinin et all d2009bh 
in calculation of Markov Chains with CosmoMC software. 
Also one can use this module to include this likelihood into 
existing parameters chains using priors adjustment proce- 
dure. In our work below we mainly used our own calcula- 
tions of Markov Chains with CosmoMC software. In some 
cases we also used chains, version WMAP7.2, taken from 
LAMBDA archive 7 . 

The constraints for some cosmological parameters in 
various cosmological models are presented in Table Q] 
More extended information on cosmological constraints 
obtained in frames of our work can be found at correspon- 
dent web page 8 . 

In all Figures below we give contours for 68% and 95% 
confidence regions. In most Figures the pair of contours of 
larger size corresponds to the set of cosmological data with- 
out the data on cluster mass function, the pair of contours 
of smaller size — with cluster mass function data taken in 
account. All confidence intervals are given at 68% confi- 
dence level, all upper limits — at 95% confidence level. 

SYSTEMATIC UNCERTAINTIES 
Systematic uncertainties of galaxy clus ter mass func 



tion me asurements were studied in detail bv lVikhlinin et al 
(l2009bl) . They were not included in likelihood functions, 



calculated in this work. Thus, systematical errors for cos- 
mological parameters constraints should be estimated sep- 
arately. These uncertainties can be estimated by measuring 
the shifts of likelihood maximum due to variation of ob- 
servables by their systematical errors. 

For galaxy cluster mass function data the main source of 
systematical errors is the uncertainty in cluster mass mea- 
surements. The mass measurements used in our work are 
based on hot intracluster gas temperatures and masses and 
were calibrated using hydrostatic measure ments of their to 



tal gra vitational masses in nearby clusters ( Vikhli nin et al 
2009ah . Systematical error of hydrostatic mass measure- 
ments was estimated as SM/M w 0.09 from the compari- 
son with the clu ster mass measur eme nts using weak lensi ng 
data taken from iHoekstral d2007l) and lZhang et all d2008l) . 

This systematical uncertainty of cluster mass measure- 
ments gives th e uncer tainty Jcrg « 0.02 for fixed fi m 



d Vikhlinin et all l2009a[) . When the other cosmological 



data are taken in account, systematical underestimation of 



http://hea.iki.rssi.ru/400d/cosm/ 
7 http://lambda.gsfc. nasa.gov/product/map/dr4/parameters.cfm 
8 http://hea. iki.rssi.ru/400d/cosm/combined/ 





Fig. 1. Constraints for Q„, and a$ in ACDM model. Contours 
of larger size — from WMAP7+BAO+Hq dataset, of smaller size — 
from CL+WMAP7+BAO+Hq dataset. Dashed lines show contours for 
SM/M rs 0.09 systematical shift in cluster mass measurements (left), and 
also for the departure from self-similar evolution, SM/M « 0.05 for dis- 
tant clusters (right). 



cluster masses for SM/M sa 0.09 leads to the shift of confi- 
dence contours, as it is shown in the left panel of Fig.Q] In 
this Figure, as an example, we show the constraints in Q. m 
- a% plane, in the model of flat Universe with cosmological 
constant (ACDM), obtained using CL+WMAP7+BAO+H 
cosmological data set. In the left part of this Figure the 
confidence contours for SM/M as 0.09 systematical shift in 
cluster mass measurements are shown with dashed lines. 
Since the likelihood maximum is shifted not along the line 
of constant f2 m , the systematic error for parameter turns 
out to be substantially smaller, Sa^ ps 0.011 in this case. 
Therefore, systematic errors for cosmological parameters 
may be significantly reduced when the additional cosmo- 
logical data are taken in account. 

In addition to the error of the mass measurement calibra- 
tion for nearby clusters, a significant part of total systemat- 
ical error comes from the uncertainties in their cosmolog- 
ical evolution, which lead to the difference in mass scale 
calibrations for nearby and distant clusters. The system- 
atic uncertainty due to the departure from self-similar evo- 
lution can be estimated as S M/M w 0.05 for distant clus- 



ters at z w 0.6 (see details in lVikhlinin et al.U2009al) . This 



difference in mass scales corresponds to the shift of con- 
fidence contours in Q, m -a% plane shown in right panel of 
Fig-Q] Quadratic sum of to systematic errors estimated for 
these two main uncertainties in cluster mass measurements 
is given below as an estimate of total systematical error for 
cosmological parameters constraints. 

In correspondent parts of this paper the systematic un- 
certainties of the data on observations of superno vae Type 



la, UNION2 compilation ( Amanullah et al. . 2010l) . are also 



considered. Fore these data the total, systematical and sta- 
tistical errors were calculated using correspondent covari- 
ance matrix. Separate estimate of systematic error for these 
data was obtained by quadratic subtraction of statistical er- 
ror from total one. 
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Table 1. Cosmological parameters constraints 



Model 


Data 


Parameter 


Value* 


Systematical error 


ACDM 


CL+ WMAP7+BAO+H 


n m h 2 


0.1311 ±0.0016 


±0.0018 








3155 ± 38 


±44 






<?% 


0.792 ±0.010 


±0.013 








0.2558 ± 0.0077 


±0.0063 






H 


71.6il.0kms" 1 Mpc" 1 


±0.4 km s" 1 Mpc" 1 








0.04386 ±0.00086 


±0.00064 


ACDM 


CL+ WMAP7+BAO+H a +SPT 


Q,„h 2 


0.1310±0.0014 


±0.0016 








3154±33 


±39 






Pi 


0.7921 ±0.0089 


±0.0115 








0.2560 ±0.0067 


±0.0072 






H 


71.55 ±0.77 km s" 1 Mpc" 1 


±0.56 km s" 1 Mpc" 1 








0.04349 ±0.00075 


±0.00071 


ACDM+Q, k 


CL+ WMAP7+BAO+H +SPT 


n* 

Qa 


-0.0018 ±0.0042 
0.7430 ±0.0072 




ACDM+m v 


CL+ WMAP7+H +BAO+SPT 




< 0.32 eV 






CL+ WMAP7+H +BAO+SPT+SN 




< 0.28 eV 




ACDM+N,,ff 


CL+ WMAP7+BAO+H +SPT 


Npff 


< 3.74 






CL+ WMAP7+BAO+H +SPT+SN 


N,ff 
e JJ 


< 3.70 




ACDM+ntv+Neff 


CL+ WMAP7+BAO+H +SPT 


Em, y 


< 0.72 eV 










n ^6 + o 16** pV 


+0 OR pV 
znu.uo c v 






Neff 


<4.62 




WCDM 


CL+ WMAP7+BAO+H 


w 


-1.027 ±0.069 


±0.028 




CL+ WMAP7+BAO+H +SN 


w 


-0.990 ±0.034 


±0.041*** 




CL+WMAP7+BAO+H +SPT 


w 


-1.013±0.066 


±0.029 




CL+ WMAP7+BAO+H Q +SPT+SN 


w 


-0.982 ±0.032 


±0.038*** 


WCDM+n k 


CL+WMAP7+BAO+H a 


w 


-1.12 ± 0.13 






CL+WMAP7+BAO+H +SN 


w 


-0.991 ±0.039 




WACDM 


CL+WMAP7+BAO+H +SN 


Wo 
W a 


-1.13 ±0.1 1 
0.47 ±0.36 





* — all intervals are given at 68% confidence level, all upper limits — at 95% confidence level; 
** — interval is non-gaussian, for significance of non-zero value see text; 
*** — including systematic uncertainty of SNe la data; 



CL+ WMAP7 CL+ WMAP7+H CL+ WMAP7+H +BAO CL+ WMAP7+H +BAO+SPT 




0.2 0.25 0.3 0.35 0.2 0.25 0.3 0.35 0.2 0.25 0.3 0.35 0.2 0.25 0.3 0.35 



^"^/JJ ^""^/ff ^""^/ff ^"^/7J 

Fig. 2. Constraints on f2„, and <rg in ACDM model, obtained using different datasets. Larger contours, from left to the right: WMAP7, WMAP7+Hq, 
WMAP7+BAO+Hq, WMAP7+BAO+Hq+SPT. Smaller contours — using the same data, with the data on galaxy cluster mass function added. 
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< 
c: 






0.125 0.13 



0.04 0.042 0.044 0.046 0.048 



Fig. 4. Constraint on spacial curvature in Universe, Q/,, in ACDM model 
with t^O. Contours are the same as in Fig. [3] 



spectral index of density perturbations n s . 

Note, that in considered cosmological model it is sug- 
gested that there are no other relativistic particles at 
equipartition in addition to photons and three known neu- 
trino species and the relativistic energy density is known 
exactly. Therefore, in this model the constraint on Vl m h 2 is 
equivalent to the constrain on equipartition redshift z eq (in 
the left panel of Fig.|3]the axis Q,,„h 2 may be changed to z eq , 
see also Table [TJ. If the number of relativistic particles is 
considered as free parameter, relativistic energy density is 
no longer exactly defined quantity and it may be measured 
from cosmological data. In this case, the data on galaxy 
cluster mass function allow to obtain significant constraints 
on the number of relativistic species (see below). 

EXTENSIONS OF ACDM MODEL 
Spacial curvature 

Constraints on spacial curvature in ACDM model with 
f2,t ^ 0, and also on the cosmological constant density are 
shown in Fig. [4] (see also Table [TJ. The constraint on 
is not improved significantly, there is only approximately 
20% improvement, as compared to the case where cluster 
data are not taken in account (see Fig.|4|. We note, that in 
this model galaxy cluster data allow to better constrain fi m 
and J7a separately. 

Primordial helium abundance 

If the data on galaxy cluster mass function are added 
to WMAP7+BAO+Ho dataset, the upper limit on primor- 
dial helium abundance, Y p , is significantly improved (see 
Fig. [5j left panel). However, the data on cluster mass 
function do not change significantly the constraints, which 
are obtained with SPT data taken in account (Fig. |5j right 
panel). From CL+WMAP7+BAO+H +SPT dataset we ob- 
tain Y p = 0.279 ± 0.028, i.e., the measurement of primordial 
helium abundance is shifted slightly to its standard value 
Y p « 0.25, which is obtained from the theoretical calcula- 
example, the constraints on the combination VL^li 1 , or on tions of primordial nucleosynthesis and from the measure- 



Fig. 3. Constraints on Q m h 2 , erg, Hq, Oj, and n s in ACDM model. Con- 
tours of larger size — from WMAP7+BAO+Hq+SPT dataset, of smaller 
size — from CL+WMAP7+BAO+H +SPT dataset. 



CONSTRAINTS FOR ACDM MODEL 

In the upper part of Table [TJ the cosmological parame- 
ters constraints for flat Universe with cosmological con- 
stant (ACDM model) are presented. The nature of main 
constraints is as follows. The constraints from galaxy clus- 
ter mass function are based mainly on its normalization and 
have their main effect on f2 m and as — the constraints 
turns out to be degenerate, these two pa rameters are re 
lated a s cr 8 = 0.813(^/0.25)"° 47 ±0.01 1 dVikhlinin et al. 



2009b). The relation between trg and CMB anisotropy am- 
plitude depends mostly on linear perturbation growth be- 
tween decoup ling epoch and z = (see, e.g., approximate 
formula from iHu. Jainl 120041) . which in turn strongly de- 
pends on O m . Therefore, the data on CMB anisotropy am- 
plitude allow to break the above tt m -ag degeneracy and 
these data analyzed jointly provide most accurate measure- 
ment of both a% and VL m . 

The relative power of various datasets to constrain 0„, 
and us parameters is shown in Fig. [2] One can see that the 
main power of constraints in ACDM model indeed origi- 
nate from the joint analysis of galaxy cluster mass func- 
tion and CMB data. The other cosmological data give only 
smaller improvement to the joint constraints on a% and Vt m . 

Parameters a% and n,„h 2 appear to be best constrained 
from joint analysis of all cosmological data considered in 
our work, which is shown in Fig. [3] The projection of these 
constraints on other parameters significantly improve their 
measurements in some cases. For example, the constraints 
are improved for parameters Ho and fib- However, the data 
on cluster mass function do not improve significantly, for 
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0.3 0.4 0.6 0.3 0.4 0. 



Fig. 5. Constraints on primordial helium abundance in ACDM model 
with free Y p . Larger contours — WMAP7+BAO+Hq dataset (left), and 
from WMAP7+BAO+Hq+SPT dataset (right), smaller contours — using 
the same data, with the data on galaxy cluster mass function added. 




-0.06 -0.04 -0.0S 0.03 



r dn s /d\nk 

Fig. 6. Constraints on tensor-to-scalar ratio and running spectral in- 
dex in correspondent extension of ACDM model. Larger contours 
— from WMAP7+BAO+H +SPT dataset, smaller contours — from 
CL+WMAP7+BAO+H +SPT dataset. 



ments of STVi 2 with standard value of effective number of 
neutrino species N e ff. 

Tensor modes and running spectral index 

It turns out, that existing cosmological data on galaxy 
cluster mass function does not improve the constraints on 
tensor-to-scalar ratio r and running spectral index of pri- 
mordial perturbations dn s /d\n(k) (see Fig.|6]l. 

CONSTRAINTS ON TOTAL NEUTRINO MASS AND 
EFFECTIVE NUMBER OF NEUTRINO SPECIES 

Massive neutrinos 

Massive neutrinos with masses of order of few 0.1 eV 
would produce significant suppression of the matter density 
fluctuation at galaxy cluster scale since they remain to be 
relativistic after equipartition and are started to be involved 
in gravitational instability growth only l ater (see, e.g., 



Doroshkevich et aU I198U iHu et all 11998b . This change 



the relation between linear perturbations amplitude erg and 
normalization of CMB anisotropy power spectrum, which 
strongly depend also on fl m . If the measurement of Hubble 
constant is added to the measurements of erg and normaliza- 
tion of CMB anisotropy power spectrum, the f2 m and total 
neutrino mass are both constrained with these data. 

Since massive neutrinos are involved in gravitational in- 
stability considerably later than dark matter particles and 
baryons, their presence could change the process of non- 
linear gravitational collapse of galaxy cluster haloes. The 
detailed analysis of this question using numerical simu- 
lations show, that these changes are indeed notable, but 
the correspondent change of cluster mass function appears 
to be not larger than few percents for neutrinos with to- 
tal mass abo ut 0.3 eV and c l usters with masses of or - 
der 10 14 M Q dBrandbvge et alll2010tlMarulri et alll201 lb . 
This is smaller that the accuracy of theoretical mass func- 
tion calculati ons used in our wo rk, which is approximately 



equal to 5% ( [Tinker et al 
The constraints on 
CL+WMAP7+H dataset 



20081) . 
total neutrino mass 
are shown in the left 



from 
panel 



of Fig. □ the upper limit is Em„ < 0.29 eV (95% c. 1., 
see also Table [TJ. In the right panel of Fig. [7] we show 
the constrains, obtained with additional available data: 
CL+WMAP7+BAO+H +SPT, in this case the upper limit 
is Em„ < 0.32 eV. If the data on SN la are added, the 
constraint is Y.m v < 0.28 eV. Therefore, the additional 
cosmological data provide almost no changes in total neu- 
trino mass constraints, as compared to CL+WMAP7+Hq 
dataset. 

From these results we see, that new cosmological data 
does not significantly improve the constraints on total neu- 
trino mass, w hich were obta i ned fro m the same galaxy 
cluster data in lVikhlinin etal.1 d2009bl) — T,m u < 0.33 eV. 
Among the other recent constraints on total neutrino mass, 
we note the upper limit £ra„ < 0.58 eV, obtained us- 
ing the data on the CMB a nisotropy spectrum an d Hub- 



ble constant measurements dKomatsu et all 1201 ll) . Also 



significant constraints were obtained from the other data 
on cluster mass function or density perturbations power 
spectrum dMalinovskv et all 120081: iMantz et all 12010b : 
Reid et allBOlOt iThomas et all l2010h . The strongest up- 
per limits on total neutrino mass was publish ed in papers , 



based on the observations of Lyman-a forest dSeliak et al 



2006b . However, the more accurate treatment of system- 
atic errors should w eaken these constraints considerably 
(Bolton et al.. 2008). A recent review of total neutrino 



m ass measuremen t s from astrophysical data can be found 



m 



Abazaiian et al.l d201 11) . 



We note, that all the constraints on total neutrino mass 
described above were obtained in assumption that only 
photons and three known neutrino species are relativistic 
at matter-radiation equipartition epoch. It turns out that if 
additional neutrino species are allowed in the model, the 
constraints on total neutrino mass change significantly (see 
below). 

Effective number of neutrino species 

The present CMB radiation density is directly related 
to the density of relativistic matter at matter-radiation 
equipartition epoch, only if the number of relativistic parti- 



ASTRONOMY LETTERS Vol.38 No. 2012 



Cosmological parameters constraints from galaxy cluster mass function 



7 




0.2 0.4 0.6 0.3 0.4 0. 



Em„, eV £m„, eV 

Fig. 7. Constraints on total neutrino mass in ACDM model with non-zero 
neutrino mass. Larger contours — from WMAP7+Hq dataset (left) and 
from WMAP7+BAO+H +SPT dataset (right), smaller contours — from 
the same data, with the data on galaxy cluster mass function added. 




Neff N e ff 

Fig. 8. Constraints on the effective number of neutrino species in ACDM 
model with free Neff. Larger contours — from WMAP7+BAO+Hq dataset 
(left) and from emphWMAP7+BAO+// +SPT (right), smaller contours — 
from the same data, with the data on galaxy cluster mass function added. 



cles species at equipartition is assumed to be known — a set 
of relativistic particles consists of photons and three known 
neutrino species. If the number of relativistic particles at 
equipartition differ from its standard value, than the rela- 
tivistic matter density is not defined anymore. In this case 
the determination of relativistic matter density at equiparti- 
tion from cosmological data gives the measurement of the 
effective number of relativistic species. 

The number of relativistic species is usually 
parametrized with the number of species of light neutrinos 
N e ff, which gives the following relation between the 
densities of relativistic mater and CMB photon energy: 



Pr = 



7 

8 V 11 



4/3 



N,„ 



Note, that for the case of three known neutrino species 
accurate calculations of neutrino decoupling give some- 
what larger effective number of neutrino species, N e ff = 
3.046 (e.g.. iGnedin & Gnedinl Il998t bolgov et all Il999 : 
Manganoetal.[l2002h . 

One of the main observables, which is measured from 
the CMB observations is the size of particle horizon at 
radiation-matter energy density equ ipartition and, there 



Hu. Dodelson. 2002 



fore, equipartition redshift z eo (e.g. 
Gorbunov & Rubakov , 2010b . On the other hand, as it was 
shown above, galaxy cluster mass function data allow to 
significantly improve the measurement of Vt„,h 2 , using the 
measurement of the matter density perturbation amplitude 
<7g and CMB anisotropy normalization. The measurement 
of two parameters, z eq and fl„,h 2 , gives the measurement of 
relativistic energy density at equipartition, which allow to 
obtain constraints on N e ff. 

The constraints on the effective number of light neutrino 
species obtained from joint analysis of data on galaxy clus- 
ter mass function and WMAP7+BAO+H^ dataset is shown 
in left panel of Fig. [8] the upper limit is Neff < 4.07. If 
the data on CMB power spectrum at higher multipoles are 
added, the constraint is improved, because the damping of 
CMB power at small angular scales depends on N e ff as 
well. This is shown in the right panel of Fig. [8] where 



SPT data are added to WMAP7+BAO+H a dataset. With 
these data taken in account, the upper limit is N e g < 3.74. 
Note, that the data on galaxy cluster mass function not only 
improve the constraints on N e g, but it is also shift them 
closer to the standard value N e g = 3.046. The data on SN 
la also improve this constraints slightly, in this case we ob- 
tain N e ff < 3.70. Systematical errors in cluster mass mea- 
surements make these upper limits slightly less strong. For 
example, if cluster masses are underestimated for SM/M = 
0.09, from CL+WMAP7+BAO+H Q +SPT dataset we obtain 
Neff < 3.89. 

These results can be compared to the following con- 
straints obtained recently. Using the results of Atacama 
Cosmology Telescope (ACT) observations in combination 
with the data on Ho and BAO measurements, the follow- 
ing constraint was o btained: N e g = 4.56 ±0.75 (68% c.l., 
Dunklev et aill2010h . Using new data of South Pole Tele- 



scope (SPT) survey together with the same measurements 
of H a nd BAO, the constra int N e g = 3.86 ± 0.42 was ob- 
tained ( Keisler et al. , 201 lb . In these constraints some- 
what higher values of N e g are preferred, but standard value 
Neff = 3.046 is not ruled out at 95% confidence level. 



In iKeisler et al.l (1201 lb . SPT data were analyzed jointly 
with the data on galaxy cluster mass function in the form 
of constrai nt on gs(^M /0.25) 47 p arameter combination, 
taken from IVikhlinin et al.l (I2009bb . In our work stronger 
constraints on N e g were obtained since we used complete 
likelihoods for gala xy cluster mas s funct ion data. Note also 
that, as compared to lKeisler et al.l (1201 lb . systematic errors 
are not included in our upper limits on N e g (see discussion 
above). 

Compatible constraints on N e g were obtained from the 
oth er data on g a laxy cl uster mass function. For example, 



in dMantz et all l2010bb the constraint N eff = 3.4^ (68% 



c.l.) was obtained. Therefore, the constraints from the 
data on galaxy clusters are in better agreement with the 
standard value N e g = 3.046, as compared to ACT and SPT 
results. A significant constraint, which is also consistent 
with the standard value N e g = 3.046, was obtained from 
the comparison of observed abundance of light elements 
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Fig. 9. Constraints on effective number of neutrino species and 
on total neutrino mass in ACDM model with free N e ff and Em„. 
Larger contours correspond to WMAP7+BAO+Hq dataset (left) and 
WMAP7+BAO+H +SPT dataset (right), smaller contours — to the same 
data, with the data on galaxy cluster mass function added. 



Fig. 10. Constraints on effective number of neutrino species and on total 
neutrino mass in ACDM model with free Neff and Sra„ . Larger contours 
correspond to WMAP7+BAO+Hq+SPT dataset, smaller contours — to 
the same data, with the data on galaxy cluster mass function added. 



with the predictions of pr imordial nucleosynthesis theory 
dMangano. SerpicoluOl lb . 

We emphasize that our constraints on N e ff, and also 
all other constraints on this parameter, discussed above 
(except the constraint from observed light element abun- 
dance), were obtained in assumption of zero total neutrino 
mass. However, if both total neutrino mass and N e g are 
considered as free parameters, the constraints on both of 
them turn to be substantially weaker. This case is consid- 
ered below. 

Effective number of massive neutrinos 

The constraints on total neutrino mass and effective 
number of neutrino species N e g in ACDM model with ar- 
bitrary number of massive neutrinos are shown in Fig. [9] 
and [10] (see also Table [JJ). Note, that if arbitrary number 
of neutrino species are allowed in the model, the upper 
limits on their total mass turn to be substantially weaker. 
And vice-versa, the assumption on non-zero total neu- 
trino mass weaken the upper limits on number of neu- 
trino species. We obtain the following upper limits from 
CL+WMAP7+BAO+H +SPT dataset: £«„ < 0.72 eV and 
N eff < 4.62. 

Therefore, the most powerful set of cosmological data 
used in our work (CL+WMAP7+BAO+H Q +SPT) does not 
exclude the existence of one additional type of neutrino. 
Moreover, for these data the maximum likelihood is shifted 
to the model with total neutrino mass about 0.4 eV and 
number of neutrino species N e jj « 4. The improvement of 
maximum likelihood for this model, as compared to the 
model with £m„ = and N eff = 3.046, is AlnL = 1.80, 
which corresponds to A^ 2 = 3.60 for 2 degrees of free- 
dom, approximately 1.4er significance. From Fig. [TOlone 
can see that the likelihood is improved mainly due to the 
assumption on non-zero neutrino mass. When massive neu- 
trinos are added to the model with free number of neu- 
trino species, the likelihood improvement corresponds to 
A\ 2 = 3.41 for one degree of freedom, which corresponds 
to approximately 1. 9a significance. 



Systematical errors of cluster mass measurements de- 
crease the significance of introduction of these new pa- 
rameters in the model. For example, if cluster masses are 
underestimated for SM/M = 0.09, then the introduction of 
massive neutrinos into the ACDM model with free N e g 
gives A\ 2 = 2.83 and significance about 1.7er. One can 
see, that systematic uncertainties in cluster mass measure- 
ments are very significant in total neutrino mass constraints 
even with existing cluster data. Therefore, in order to im- 
prove similar measurements in future one will need to sig- 
nificantly reduce systematic uncertainties in cluster mass 
measurements. 

Different mass distributions of neutrino species 

All constraints on total neutrino mass, discussed above, 
were obtained in assumption that three known neutrino 
species have the same equal masses (and additional neutri- 
nos are massless), i.e. the mass distribution of these three 
known neutrino species is degenerate. However, from the 
observations of neutrino oscillation it is known, that masses 
of different neutrino species are also different. Generally, if 
the mass distribution of neutrino species is non-degenerate, 
the constraints on total neutrino mass from cosmologi- 
cal data may ch a nge. This happens for several reasons 
(see, e.g., Slosar , 2006t) . in particular since most massive 
neutrino species become non-relativistic earlier and spec- 
trum of linear density perturbations is suppressed at smaller 
scales. Also, if neutrino mass is larger then approximately 
0.5 eV, they become non-relativistic before recombination 
and significantly change angular power spectrum of CMB 
anisotropy. 

In case of N e g = 3.046, the assumption on non- 
degenerate neutrino mass distribution makes almost no 
change in the neutrino mass constraints from cosmo- 
logical data. For example, if there is only one mas- 
sive neutrino, the constraint on total neutrino mass from 
CL+WMAP7+BAO+H +SPT dataset weaken only very 
slightly, SbIj, < 0.33 eV (as compared to £m„ < 0.32 eV, 
in case of three degenerate neutrinos, see above). This is 
quite consistent with what is expected for similar cosmo- 
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line neutrino oscillations searches are better explained if 
two additional light sterile neutrino specie s are suggested 



Neff 



Fig. 11. The same constraints as in Fig. [5] but in assumption that there is 
only one type of neutrino with non-zero mass. 



logical data (e.g jTakada et all I2006L ISlosarl l2006h . 

It turns out, that with increasing number of neutrino 
species, the assumption that only one neutrino is massive 
weaken the constraint on total neutrino mass more signif- 
icantly. This is shown in Fig. QT| — the upper limit for 
CL+WMAP7+BAO+H +SPT dataset is Y,m v < 0.88 eV in 
this case. The upper limit on number of neutrino species 
changes only slightly and is N e g < 4.68. These upper lim- 
its can be compared to Em^ < 0.72 eV and N e jj < 4.62 in 
case of three degenerate neutrinos (see above). 

Light sterile neutrinos 

The constraints on Em„ and N e g discussed above 
can have important application to study the possibil- 
ity of existence of light sterile neutrinos with masses 
near 1 eV, which were suggested in order to ex- 
plain the results of short bas eline neutrino oscilla tions 
searches ex periments LSDN ( Aguilar et all 2001 ) and 
MiniBooNE dAguilar-Arevalo A. A. et al. . 2010h . and also 
recently discovered, so-called reactor neu trino anomaly 



(Mue ller et all 1201 lb Mention et all 1201 lb . In order to 



explain these results a possibility of th e existence of one 



or two light sterile neutrino species (Maltoni, Schwetz 


20071 Karagiorgi et a 


1 20091 Melchion 


-i et all 2009 


Akhmedov. Schwetzl 


20 10t Koppetal.. 


201 lh. which 



should be thermalized in early Universe (see, e.g., dis 



cussion and references in Hamann et al. , 2010l) was dis 



cussed. The presence of these additional neutrinos should 
also change the value of Neff, measured from cosmological 
data. 

The upper limit N e /f < 3.70, which was obtained in our 
work in assumption on zero total neutrino mass (see above), 
does not agree with the existence of even one additional 
type of neutrino. However, in order to explain the re- 
sults of short baseline neutrino oscillations experiments ad- 
ditional species of massive neutrinos a re suggested (e.g., 



Mention et all 1201 It iKopp et all 1201 II) . In this case the 



constraints on N e g turn to be significantly weaker. As it was 
shown above, all considered cosmological data are consis- 
tent with the existence of one additional type of light sterile 
neutrino with mass about 0.4 eV. 

It was discussed recently, that the results of short base- 



( Ko pp et all 1201 It iGiunti. Lavederl 1201 lal) . However, re 



cent improvements of the MiniBooNE experiment data, 
probably will allow to explain all the data of these ex- 
periments with only on e additional light sterile neutrino 
dGiunti. Laveden, 1201 lbl) . Anyway, in order to explain the 
results of these experiment assuming only one additional 
type of neutrino, the masses near 1 eV are proba bly re- 
quired ( Kopp et al. . 201 1 : Giunti. Lavedei . 201 lbl) . This 
value disagree with upper limits on total neutrino mass, ob- 
tained above — even in most conservative case of one ad- 
ditional massive neutrino the upper limit on total neutrino 
mass is Y,m u < 0.88 eV at 95% confidence level. There- 
fore, we conclude that masses of light sterile neutrinos, 
which may explain current results of short baselines neu- 
trino oscillation searches are not in good agreement with 
our constraints. 

In order to reconcile neutrino masses near 1 eV with 
the constraints on total neutrino mass from cosmological 
data it was suggested to consider the cosmological mod- 
els wi th dark energy more general than cosmological c on- 
stant dKristiansen. ElgaroviuOl ltlHamann et all 1201 lb . It 



was found that with dark energy equation of state param- 
eter w < -1 the limits on neutrino mass are relaxed con- 
siderably. This remains true also when the galaxy cluster 
mass function cosmological data taken in account. The up- 
per limits on total neutrino mass and effective number of 
neutrino species in the model with free w, and free number 
of massive neutrinos are Em„ < 0.85 eV and Neff < 4.33 in 
assumption that three neutrino species have equal non-zero 
mass, and £ra„ < 1.00 eV and N e ff < 4.39 in case if there 
is only one type of massive neutrino in the model. There- 
fore, the limits on total neutrino mass, obtained with galaxy 
cluster mass function data taken in account, are still in poor 
agreement with suggested sterile neutrino mass near 1 eV, 
even in cosmological models with dark energy equation of 
state parameter w < — 1. 

Therefore, all available cosmological data are consistent 
with one additional light sterile neutrino. However, their 
total mass may be near 0.4 eV, while the values of total 
neutrino mass near 1 eV are in poor agreement with exis- 
tent cosmological data. We emphasize that the constraints 
obtained in our work apply only to light sterile neutrinos 
thermalized in early Universe. 

DARK ENERGY CONSTRAINTS 
In order to obtain most powerful constraints on dark en- 
ergy equation of state, the measurements of accelerated ex- 
pansion of Universe from th e observations of dista nt super- 
novae type la were used in iKomatsu et al.l (1201 lh . On the 
other hand, during the last years it become clear that the er- 
rors of cosmological measurements bas ed on SN la data are 
dominated by system atic uncertainties dHicken et al. i l2009t 
Kessler et al. . 2009). In recent works on SN la observa- 
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Fig. 12. The constraints on dark energy density and equation of state 
parameter in WCDM model. Larger contours — for WMAP7+BAO+Hq 
dataset, smaller contours — for WMAP7+BAO+Hq+CL dataset (left 
panel) and for WMAP7+BAO+H +SN dataset (right panel). Dashed line 
in left panel show the contours for systematic SM/M Ri 0.09 shift in clus- 
ter mass measurements. In the right panel the contours for SN la data with 
systematic uncertainties taken in account are shown with dashed line. 
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tions, systematical errors are estimated to be very signif- 
icant and to be larger than statistical errors in these data. 
These uncertainties are mostly consist of uncertainties of 
photometrical calibrations, color corrections on the extinc- 
tion in host galaxies, selection effec ts in SN la surveys and 
others (e.g., Amanullah et al. . 2010h . 

With these considerations in mind, it is interesting to 
study the possibilities to obtain powerful dark energy con- 
straints independently on SN la measurements. We show 
below that these constraints may be obtained using galaxy 
cluster mass function measurements, and their uncertain- 
ties are not currently dominated by systematics. The up- 
dated dark energy constraints obtained from the combina- 
tion of recent cosmological data are also given below. 

Flat Universe 

The constraints on the density and equation of state pa- 
rameter w of dark energy in a fiat Universe with free w 
(WCDM model) are shown in Fig. [12] and [13] (see also 
Table Q]). With no systematic uncertainties taken in ac- 
count, the existing data on galaxy cluster mass function 
give somewhat weaker constraints, as compared to the 
data on SN la. However, the error on w parameter mea- 
sured from from SN la data is dominated by systematic 
uncertainties (see right panel of Fig. [12}. With systemat- 
ical errors taken in account, the constraint on dark energy 
equation of state parameter from WMAP7+BAO+Hq+SN 
dataset is w = -1.068 ± 0.077. On the other hand, the 
WMAP7+BAO+H +CL dataset gives the following con- 
straint: w = -1.026 ±0.069 (stat.) ±0.028 (sys.), i.e., the 
total error is ±0.074 in this case. Therefore, the constraints 
from galaxy cluster mass function data and from SN la ob- 
servations are comparably powerful. Note that, in contrast 
to SN la data, the errors from cluster mass function mea- 
surements are not currently dominated by systematic un- 
certainties. 

The constraints for all the data combined are shown in 
Fig. Qj] One can see that galaxy cluster mass function 



Fig. 13. The constraints on dark energy density and equation of state pa- 
rameter in WCDM model. Larger contours — for WMAP7+BAO+Hq+SN 
dataset (dotted lines) and also for WMAP7+BAO+H +CL (dashed lines), 
smaller contours — for all the data (WMAP7+BAO+H +SN+CL) com- 
bined. 





Fig. 14. The constraints on the curvature of Universe and dark energy 
equation of state parameter in WCDM model with ^0. Larger contours 
— for WMAP7+BAO+Hq dataset, smaller contours — for the same data 
with the data on clusters (left panel) and SN la (right panel) added. The 
contours for SN la data with systematic uncertainties taken in account are 
shown with dashed lines in right panel of the Figure. 



data significantly improve the dark energy equation of state 
parameter constraints. The reason is that these data are 
independent and have different degeneracies in parameter 
space. From these data we obtain the following measure- 
ment: w = -0.990 ±0.034 (stat.) ±0.041 (sys.), where SN 
la systematic uncertainties are also included in resulting 
systematic error. These cons traints appears to be som ewhat 
better than those obtained in lVikhlinin et alj d2009bl) . 



Curved Universe 

The constraint on the curvature of Universe and dark 
energy equation of state parameter in WCDM model with 
Clk y are shown in Fig. [14] and [15] In model with non- 
zero space curvature the cluster data allow to obtain only 
less powerful constraint on w, as compared to supernovae 
data (cf. left and right panels in Fig.fPfl). However, taking 
in account systematic errors of supernovae data, the accu- 
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Fig. 15. The constraint on the curvature of Universe and dark en- 
ergy equation of state parameter in WCDM model with f2j ^/ 0. 
Larger contours — for WMAP7+BAO+H +SN (dotted lines) and for 
WMAP7+BAO+Hq+CL dataset (dashed lines), smaller contours — for 
all the data (WMAP7+BAO+H +SN+CL) combined. 



racy of w measurement appears to be comparable. From 
Fig. Q3] one can see that the data on galaxy cluster mass 
function combined with all other data, including the data 
on SN la, allow to significantly improve the measurement 
of the parameter of dark energy equation of state. 

Time -dependent equation of state 

The constraints on the parameters of time-dependent 
dark energy equation of state, defined as: 

w(a) = wq + w fl (l — a) , 

are shown in Fig. [16] (see also Table [T] WACDM model). 
In order to calculate CMB anisotropy power spectra we 



used Parametrized Post-Friedmann approach dFang et al 
2008b . which was done using the correspondent module 
for CosmoMC software. One can see, that galaxy clus- 
ter mass function data significantly improve the overall 
constraints for wq and w a parameters. Therefore, existent 
cluster mass function data give comparably powerful con- 
straints for these parameters, as compared to supernovae 
data. 

SUMMARY 

In our work we present the results of detailed analysis of 
cosmological parameter constraints which were obtained 
from combination of galaxy cluster mass function measure- 
ments (Vikhlinin et al., 2009a,b) with the other cosmolog- 
ical data, obtained recently. We show, that this allow to 
significantly improve the constraints for many cosmologi- 
cal parameters. 

All considered cosmological data combined together are 
consistent with the model of flat Universe with cosmolog- 



Fig. 16. The constraints on dark energy equation of state parameters in 
WACDM model. Larger contours — for WMAP7+BAO+H +SN dataset, 
smaller contours — for WMAP7+BAO+H +SN+CL dataset. 



ical constant (ACDM). In frames of this model all consid- 
ered cosmological data, with galaxy cluster mass function 
data taken in account, give the most powerful constraints 
on 08 and fl m h 2 parameters. The constraints on other pa- 
rameters, such as O m , Hq, ri;,, are also improved. The data 
on galaxy cluster mass function allow to constrain these pa- 
rameters with as high as k 1% accuracy (see Table [TJ. At 
the same time, the systematical errors from the uncertain- 
ties of cluster mass measurements appear to be comparable 
to statistical ones. 

The data on galaxy cluster mass function allow to ob- 
tain new important constraints on total neutrino mass Y>m v 
and effective number of neutrino species N e ff. When mas- 
sive neutrinos or additional number of neutrino species 
are added to ACDM model, the strongest constraints are 
obtained: Em„ < 0.32 eV and N e f < 3.74. However, if 
both total neutrino mass and effective number of neutrino 
species are considered as free parameters in ACDM model, 
i.e., the model with arbitrary number of massive neutrinos 
is considered, the constraints turn out to be significantly 
less strong: T,m u < 0.72 eV and N e g < 4.62. Moreover, 
all considered cosmological data are somewhat better fitted 
by the model with non-zero neutrino mass Em„ w 0.4 eV 
and larger than standard value of the number of neutrino 
species, N e g f=s 4. 

These constraints can have an important application to 
study of the possibility of the existence of light sterile 
neutrinos with masses near 1 eV, which were suggested to 
explain the results of short baseline neutrino oscillations 
searches experiments LSDN and MiniBooNE, and also 
recently discovered, so-called reactor neutrino anomaly. 
However, as compared to cosmological constraints ob- 
tained in our work, in order to explain the results of 
these experiments, too large number of additional neu- 
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trino species (more than one) or too large neutrino masses 
Eraj, ps 1 eV are required. 

In our work the updated constraints on dark energy equa- 
tion of state parameters are also presented. The constraints 
obtained using galaxy cluster mass function data, with no 
use of SN la measurements, are considered as well. We 
show, that taking in account systematical uncertainties, 
dark energy parameters constraints from cluster mass func- 
tion data are comparable in their power with those from SN 
la observations. Moreover, in contrast to SN la data, the er- 
rors of dark energy parameters measurements from current 
cluster mass function data are mainly statistical, and are not 
dominated by systematic uncertainties. Therefore, expan- 
sion of a sample of galaxy clusters, suitable for accurate 
measurements of galaxy cluster mass function, will allow 
to further improve dark energy constraints in future. 
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